Coherent transport in graphene nanoconstrictions 
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We study the effect of a structural nanoconstriction on the coherent transport properties of 
otherwise ideal zig-zag-edged infinitely long graphene ribbons. The electronic structure is calculated 
with the standard one-orbital tight-binding model and the linear conductance is obtained using the 
Landauer formula. We find that, since the zero-bias current is carried in the bulk of the ribbon, this 
is very robust with respect to a variety of constriction geometries and edge defects. In contrast, the 
curve of zero-bias conductance versus gate voltage departs from the (2n + l)e 2 /h staircase of the 
ideal case as soon as a single atom is removed from the sample. We also find that wedge-shaped 
constrictions can present non-conducting states fully localized in the constriction close to the Fermi 
energy. The interest of these localized states in regards the formation of quantum dots in graphene 
is discussed. 
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I. INTRODUCTION 



The recent fabrication of field effect transistor devices 
based both upon quasi-2D graphite quantum dots 1 - and 
upon graphen o 2 i 3 i 4 (a single atomic layer of graphite), 
and the observation of a new type of quantum Hall ef- 
fect in the latter have triggered a huge interest in the 
electronic properties of this system. Most of the results 
of standard mesoscopic physics need to be revisitec&Sii 
in the case of graphene because its electronic structure is 
fundamentally different from that of metals and semicon- 
ductors where either a large density of states (DOS) at 
the Fermi energy or a large gap determine the properties 
of the materials. Graphene is a semimetal with zero DOS 
at the Fermi energy and zero gap. On top of that, the 
electronic structure close to the Fermi energy has a con- 
ical shape with perfect electron-hole symmetry, identical 
to that of two dimensional massless Dirac fermions 8 . 

Here we consider graphene-based one dimensional flat 
structures, the so called graphene nano-ribbons&ifliii* 1 ^. 
As opposed to two dimensional graphene, ribbons can 
present a finite density of states at the Fermi en- 
ergy which dominates electrical transport in undoped or 
weakly doped samples. Ideal graphene ribbons can be 
considered as the flat parent structures of carbon nan- 
otubes, whose electronic properties have been thoroughly 
studied- 4 ". Electronic transport in carbon nanotubes has 
been studied in different regimes, including ballistic^, 
Coulomb Blockade 16 and Kondoiii 8 .. Pro gress in the 
fabrication of graphene based nano-structures that per- 
mits to study transport in graphene ribbons motivates 
this work. 

The presence of edges makes the electronic structure of 
graphene ribbons different from that of nanotubes. Two 
types of idealized edges are usually considered: Arm- 
chair and zig-zag 9 . Interestingly, all the zig-zag and 
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some of the armchair edges result in a band at the Fermi 
energy^. In the case of narrow zig-zag-edged ribbons, 
the top and bottom edge states can be sufficiently close 
as to make hybridization possible, resulting in two low 
energy dispersive bands, symmetrically placed around 
the Fermi energy, E = . Depending on their width, 
nano-ribbons with armchair edges can be metallic or in- 
sulating. The different density of states for armchair and 
zig-zag edges has been experimentally observed in Scan- 
ning Tunnel Microscope (STM) experiments with atomic 
resolutioni2i2S. 

Coherent or quantum transport in graphene rib- 
bons has been studied previously, both for ideal 1 ^ and 
defectiv o 21 i 22 cases in the infinite length L case. The 
case of disorder free graphene ribbons with finite L com- 
parable to the width W has also been studied in con- 
nection with the experimental measurement of a mini- 
mum conductivity in mesoscopic size graphene layers^. 
In the spirit of quantum point contact physics, in this 
work we study the coherent transport of infinitely long 
graphene ribbons, narrower than W ~ 5nm, with a struc- 
tural nanometric constriction like the one schematically 
shown in Fig. (JTJ . This type of structure can also be con- 
sidered as an idealization of an all-carbon single-molecule 
junction. In conventional single-molecule junctions the 
organic molecule is attached to metallic electrodes. Here 
these are replaced by perfect graphene ribbons and the 
role of the molecule is played by a geometrical constric- 
tion. Due to the very different electronic structure dis- 
played by different all-carbon nano-structures 1 ^ the con- 
duction properties of these systems do not appear obvious 
a priori to us. 

We calculate the electronic structure in the one orbital 
tight binding approximation. The relevant orbital is the 
p z , since the sp 2 orbitals form bonding and anti-bonding 
states very far away in energy. In the TB approximation 
ideal two dimensional graphene has a conical energy dis- 
persion at low energies and the Fermi surface (E = 0) 
is formed by six points. This low energy region can be 
described in terms of a k ■ p theory whose mathematical 
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structure is very similar to the Dirac theory for massless 
fermions 8 , which yields physical insight. The k ■ p has 
also been worked out for edge states in ideal ribbons^. 
The TB approach provides natural energy and momen- 
tum cutoffs to the k ■ p theory permits to model per- 
turbations at the atomic scale and is a good preliminary 
step towards ab-initio calculations^. 

The rest of the paper is organized as follows. In section 
II we review the transport formalism. In section III we 
review the electronic structure of ideal ribbons and we 
study the effect of vacancies, located at the edge of the 
ribbons, on their transport properties. We find that the 
conductivity of undoped graphene ribbons is weakly af- 
fected by the presence of this kind of disorder. In section 
IV and V we study transport properties of square- and 
wedge-shaped constrictions. Whereas the former present 
finite conductance, the latter have a vanishing transmis- 
sion at low energies, coexisting with a finite density of 
localized states. These states at energies close to zero 
form what can be called a quantum dot in graphene. In 
section VI we discuss the validity of our approximations 
and we summarize and discuss the main results of this 
manuscript. 

II. TRANSPORT THEORY 

The natural framework for transport calculations in 
nanoscopic devices is the Landauer formalism. The de- 
scription of electron transport within the Landauer for- 
malism is based on the assumption that transport across 
the highest resistance region is coherent, i.e., inelastic 
scattering is negligible there. A more complete account 
can be found, e.g., in the book by Datta2i. In what fol- 
lows we will assume that the inelastic mean free path for 
graphene electrons is much longer than the typical di- 
mensions of the nanoconstrictions considered and than 
the ribbon widths. According to recent theory work, in- 
elastic scattering due to phonons is very inefficient in 
graphene ribbons so that our assumption seems to be met 
even at room temperature^. This makes elastic scatter- 
ing the major contributor to resistance. In contrast to 
metallic systems where the resistance of the electrodes 
is negligible compared to that of the nanoconstriction, 
in the case of graphene nanoconstrictions, the low- bias 
conductance of the ideal electrode is also very small, =f-. 
The consequences of this are explored in what follows. 

We consider the effect of constrictions on the trans- 
port properties of an otherwise ideal ribbon (Fig. 
The system has three regions: the central region (or de- 
vice) where the constriction is located and the left and 
right leads. The latter are described as semi-infinite one- 
dimensional perfect ribbons of finite width, characterized 
by the number of atoms in the unit cell, N. We only con- 
sider ribbons with zig-zag edges as the ones studied in 
previous worka&ifiiSi. We consider constrictions of vari- 
ous shapes, from the removal of a single atom or a few 
atoms on the edge to square-shaped and wedge-shaped 
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FIG. 1: Division of system into leads L and R, device D, and 
scattering region S (LEFT) . And division of operator matrices 
into corresponding submatrices (RIGHT). 

constrictions. 

The introduction of the constriction breaks the trans- 
lational symmetry of the perfect lead so that, in general, 
electrons incident on a given band will be either reflected 
or transmitted into other bands after hitting the constric- 
tion. The square of the transmission amplitude t nm gives 
the probability of an incoming mode m to be scattered 
on an outcoming mode n. The Landauer formula links 
the overall transmission 

T{E) = IIW|| 2 - Tv[th] 

with the linear conductance 

G{E F ) = ^-T(E F ) 
h 

where Ep is the Fermi energy. 

For completeness, we review the basic steps in the cal- 
culation of the transmission T(E) using single-particle 
Green's functions as routinely done in nanoelectronics^. 
In the TB approach used here, the Hamiltonian matrix is 
straightforwardly obtained for a given atomic structure. 
The leads are characterized by a unit cell with N atoms 
and a propagation direction which we take along the x 
axis. When written in blocks N x N, the Hamiltonian of 
the leads is a semi-infinite tridiagonal matrix, with intra- 
cell blocks Hq and an inter-cell first-neighbor coupling 
V. Our starting point is the partition of the Hamilto- 
nian of the infinite system in 3 regions. The choice of the 
boundaries between device and leads is done so that left 
and right electrodes are not directly coupled. For con- 
venience the device region is always chosen so that the 
left and right boundaries are given by one of these units 
cells. Therefore, the coupling between the surface unit 
cell of the left (right) electrode to the right (left) bound- 
ary of the device is given by the same inter-cell matrix V. 
According to this scheme, the Hamiltonian matrix set is 
divided into submatrices as follows: 

/ H L H LB \ 
H = [ H nL H D H DR , (1) 
\0 H RD H R ) 

Here Hy> is a finite size square matrix with range equal 
to the number of atoms in the device, N4. In contrast, 



Hl.r are infinite size square matrices describing the semi- 
infinite electrodes and i?DL and i?LD are infinite rectan- 
gular matrices. The Green function operator, defined as 



(E - H)G(E) = 1 
can also be divided into submatrices as: 

(Gl Gld Glr 
Gdl Gd Gdr 



(2) 



(3) 



After simple steps, it is possible to write the Green 
function of the device as: 



Gu(E) = (E-H B - E L (£) - £ R (.E))- 



(4) 



where El,r are the so called self-energy Nd X Nd matrices 
given by 



£l(£) := H DL g h {E)H hD 
E K (E) :=H DR g R (E)H R u 



(5) 



The selfenergies Y,l(E) and Er(E') describe the effect 
of the electrodes on the electronic structure of the de- 
vice. The real part of the self-energy results in a shift 
of the device levels whereas the imaginary part pro- 
vides a lifetime. The device self-energies are given by 
the Green's functions of the semi-infinite isolated leads 
g L (E) = (E-Hl)- 1 &ndg K (E) = (E — Hr)^ 1 projected 
into the device region by the coupling of the leads to the 
device 7?dl and Hrd- In contrast to the Green func- 
tion of an infinite system with translational invariance, 
the calculation of the Green function of a semi-infinite 
system with a surface is non-trivial. 

We can write the surface part of the semi-infinite Green 
function <7l,r(^) of the electrode as: 



Sl,r(£)| s 



1 



rfacc 



E — Hq — Ei jr (E) 



where Ei >r is a self-energy (different from El,_r) that ac- 
counts for the effect of the rest of the semi-infinite chain 
on the first unit cell. In one dimension it is possible to 
derive a recursive relation that yields a self-consistent 
Dyson equation for this self-energy: 



El,rCE) = V 



1 



E — Hq — Ei, r (£7) 



The coupling matrices Ti,(E) and T R (E) are defined as 
the difference between the retarded and advanced selfen- 
ergy of the leads projected into the device by the coupling 
if DL and -Hrd : 



T(#)l(r) = i ( S L(R) 



(E) 



J L(R) 



(E) 



With all these ingredients, we can compute the trans- 
mission using the result^: 

T(E) = Tr [Tl (E)G\^ (E)T r (E)Gy> (E)] . 
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FIG. 2: (a) Bands of the ideal zig-zag ribbon with 24 atom 
unit cell, (b) Detail of the low energy sector, (c) Detail of the 
low energy sector for the 48 atom unit cell 



Therefore, for a given system, we first compute the elec- 
trode surface Green function © by solving the Dyson 
equation (7J|. This permits to compute the device self- 
energy (0, the device Green function and the cou- 
pling matrices ((HJ- The final step is the calculation of 
the transmission function. 



III. IDEAL AND WEAKLY DEFECTIVE 
ZIG-ZAG RIBBONS 

A. Ideal ribbons 



In this section we briefly review the electronic structure 
and transport properties of both ideal and weakly defec- 

(6) tive zig-zag-edged ribbons. We consider ribbons with two 
different widths whose unit cells are composed of N = 24 
and N = 48 atoms, respectively. The super-cell unit of 
the N = 24 ribbon is shown in Fig. (J3h>) ■ The honeycomb 
lattice of 2D infinite graphene can be generated by a tri- 
angular lattice of unit cells with two atoms, labeled A and 
B. Therefore, the honeycomb lattice is formed by two in- 
terpenetrating sub-lattices, A and B. The first neighbors 
of atoms in the lattice A belong to sub-lattice B and vice 

(7) versa. This underlying structure is responsible for most 
of the peculiar features of graphene electronic structure. 
Nano-ribbons inherit these properties and, in the case of 
zig-zag nano-ribbons, the top and bottom edges belong 
to atoms on different sublatticesi^. 

Figure (0) shows the band structure for the N — 24 
case [Figs. |2Ja) and^fb)] and for the N = 48 case [only 

(8) the low energy region is now shown in Fig. Hfc)]. The 
energy units are given in terms of the hopping parameter 
\t\ ~ 3eV which is the only energy scale in the Hamil- 
tonian. There is a perfect electron-hole symmetry which 
sets the Fermi energy to zero for half filling. Notice that 

(9) the density of bands per energy interval is largest for the 
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FIG. 3: Upper panel: Local density of states of the ideal 24 
atoms-wide ribbons for E — (circles) and E = O.Oli (tri- 
angles). Notice the vertical logarithmic scale. Lower panel: 
atomic structure of a a section of a N = 24 zig-zag ribbon. 



wider ribbon, as expected. We can distinguish three dif- 
ferent regions, (i) \E\ < Ai, (ii) Ai < \E\ < \t\ and (hi) 
\E\ > t, where Ai is the minimum energy of the second 
sub-band closest to \E\ = 0. This energy scale Ai is asso- 
ciated to the finite width of the ribbon and it decreases as 
N^ 1 . For the cases considered here, N = 24 and N = 48, 
Ai = 0.35 t and Ai = 0.18t respectively. In the N = oo 
limit Ai goes to zero, as expected for two-dimensional 
graphene. 

The density of extra carriers that needs to be injected 
in the ribbon so that the Fermi energy hits the second 
subband, Ep = Ai, exceeds the upper experimental limit 
reached so fa»i, Sn — 10 13 cm -2 . Electrical doping up to 
the second subband could be possible for wider ribbons. 
Although only the lowest energy region \E\ < Ai could 
be accessible experimentally for the narrow ribbons con- 
sidered here, some insight can gained also by analyzing 
the effect of elastic scattering on the transmission at en- 
ergies in the other two regions. The bands immediately 
above Ai present two positive (and two negative) mo- 
menta at a given energy. This results in the doubling of 
the number of conducting channels at a given energy, as 
long as \E\ < t. Higher in energy, \E\ > t, we find simple 
parabolic bands which yield one channel per band as in 
the case of III-V semiconductors. 

The energy bands closest to Ep are flat for \k\ larger 
than a critical wave- vector k Cl and present dispersion oth- 
erwise. The cutoff k c is a decreasing function of N which 
tends to k c a/n — 2/3 in the N — oo limit 9 . In figure[31we 
plot the local density of states across a super-cell unit for 
the ribbon with N — 24 and for two different energies, 
E = and E = O.Olt. Atoms A and B in the same unit 
cell are joined by a straight line. In both cases the LDOS 



is peaked on the edges, although the contrast is larger for 
the E = case, which corresponds to the dispersionless 
states. The LDOS presents a peculiar oscillation, related 
to the sub-lattice structure. Starting from the left edge, 
the LDOS peaks only in A, the weight on the first B 
atom being very small. As we move towards the center 
the weight on A atoms decays exponentially whereas the 
weight on the B atoms increases exponentially. In the 
center the weight on the A and B atoms is comparable 
and the weight on the B atoms becomes dominant as 
we move towards the opposite edge, in very good agree- 
ment with the k ■ p theorjii^. As a result, the density 
of states, disregarding the sub-lattice index, is peaked at 
the edges, which permits to refer to the lowest subband 
states as " edge states" . 

It is apparent from Fig. [3]that the A and B edge states 
become coupled in the middle of the ribbon. Importantly, 
the structure of the current operator is the same as the 
hopping part of the Hamiltonian, coupling atoms of the 
different sub-lattices. Therefore, the current density of a 
given state, evaluated in a unit cell with two atoms, is 
proportional to the product of its A and B components, 
ipA{y) and iI>b(jj)- hi the continuum limit the current 
associated to the wave function ^ — (V^jV'b) rea ds 
jx(y) oc C&<cr x ty), where a x is the Pauli matrix acting 
on the AB spac e 5 ' 6 ' 26 . From figure (J3J) we expect that 
the product ipA{y) x ^Pb{v) to take similar values in the 
edges 12 and in the center of the ribbon. The resulting 
picture is the following: whereas the charge density of 
low energy states, proportional to p(y) — \iPa\ 2 + \iPb\ 2 , 
is peaked on the edges, their current density j x (y) is more 
homogeneously distributed accross the ribbon. This pic- 
ture is substantially different from that of non-relativistic 
electrons with scalar wave functions <f>{x,y) = il}(y)e lkx 
for which the charge density p{y) = |t/>(?/)| 2 and the cur- 
rent density j x (y) oc <pd x 4>* — <p*d x <j) oc kp(y) have the 
same profile. The consequences of this difference between 
the current and charge densities will become apparent 
later and are one of the results of this work. 



B. Defective edges 

We now study the effect of a single vacancy in the 
edge(s) of the ribbon on the transport properties. From 
the formal point of view this is done using the approach 
described in the previous section. The sector of the rib- 
bon where the vacancy is located is treated as the device. 
In Fig. 2| we plot T{E) for three cases: Ideal case, one 
atom missing and two atoms missing. As everywhere 
else in the text, this is the transmission per spin chan- 
nel. Since we assume that time-reversal symmetry is not 
broken, the total transmission should be multiplied by a 
factor of 2 to account for the spin degree of freedom. All 
of them display electron- hole symmetry around E = 0. 
As in the case of the energy bands, the T(E) curves have 
three different regions. In the large \E\ region, the trans- 
mission for the ideal ribbon T l( icsd(E) is quantized ac- 
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FIG. 4: (Color on-line), (a) Transmission spectrum for rib- 
bons with N = 24 atoms in the ideal case (black) , removing 1 
atom from the edge (red), and removing 2 atoms on the top 
and bottom edges (blue), (b) Detail of the transmission in 
the low energy sector. The effect of removing I or 2 atoms on 
the border is very small in the lowest conductance plateau. 



cording to the usual law for non-relativistic electrons and 
holes: ne 2 /h, where n is the number of bands that inter- 
sect the Fermi energy for k > 0. In this case n goes from 
1 to 12 in steps of 1, consistent with the parabolic dis- 
persion away from the Dirac cones. This large \E\ region 
is not likely to be reached experimentally since it would 
imply a huge depletion of the charge, but is interesting 
from the conceptual point of view. In the intermediate 
region (t < \E\ < Ai) T- lt \ ca \(E) changes according to 
the rule (2n + l)e 2 /h. This result has been previously 
obtained by Peres et al. lfl3i) . The factor of two comes 
from the shape of the bands in this region, as discussed 
above. The transmission is maximal at the energy where 
the "Dirac" ladder and the non-relativistic ladder meet 
at E ~ /. In the relevant low energy sector, where we 
will focus our attention from now on, the transmission of 
the ideal ribbon is one. 

The effect of scattering introduced in the ribbons af- 
fects very differently the transmission in the different en- 
ergy regions. The removal of a single atom on the edge 
erases the 2n + 1 ladder in the intermediate energy re- 
gion for the N — 24 ribbon [see Fig. Ufa)]. This points 
towards a very difficult experimental verification of the 
In + 1 transmission ladder in these narrow ribbons. In 
contrast, in the case of N = 48 the conductance remains 
quantized in the 2n + 1 = 3 plateau after the removal 
of 1 atom in each edge, but the higher energy plateaus 
disappear. 

In contrast with the higher energy subbands (\E\ > 
Ai), the plateau in the low energy sector is only weakly 
affected by the removal of atoms in the edges, even for 
the narrower ribbon (N = 24) as can be seen in Fig. 
IHb). The effect on the transmission in the experimen- 
tally relevant low energy sector is below two percent for 
the removal of one atom on one edge and below five per- 



cent when one atom is missing on both edges. In contrast, 
the removal of a single atom in the central part of the rib- 
bon (not shown) has a much larger influence on the low 
energy transmission. These results are compatible with 
the fact that the current density carried by edge states 
is spread along the central region. The robustness of the 
transmission in the low energy sector and its weakness in 
the intermediate energy sector are relevant results since 
atomic size fluctuations in the edges are unavoidable in 
real samples. 



IV. SQUARE CONSTRICTIONS 

We now present results for square-shaped nanocon- 
strictions as those shown in Fig. |SJ We choose this shape 
because it permits a comparison with a simple model for 
electrons in parabolic bands. It also seems possible to ob- 
tain analytical expressions for the transmission curve for 
square constrictions using analytical approach along the 
lines of previous work^iiS. Another feature of square 
constrictions is the presence of armchair edges joining 
the zig-zag edges of the wide and narrow regions. We 
consider square constrictions with top-bottom symme- 
try, like those in the Fig. (JSJ, characterized by three 
lengths: The width of the electrode and the width of the 
constriction, N and N c , both measured in units of the 
number of atoms of the super-cell unit, and the length of 
the constriction L, measured in units of a, the graphene 
lattice parameter. 




FIG. 5: (Color On-line). Structure of the square-shape 
nanoconstriction with 3 characteristic lengths, N, N c and L. 

In Fig. @ we show transmission curves for a variety 
of square constrictions. In Figs. Efa) and (b) we show 
results for electrodes with N = 24 atoms, whereas in^c) 
we show results for wider electrodes with N — 48. The 
two curves in Fig. EJa) permit to compare the trans- 
mission for two structures with the same aspect ratio 
N c /N = 12/24 but different nanoconstriction length L. 



Analogously, Fig. |(Jb) shows two T(E) curves with the 
same aspect ratio N c /N — 8/24 and different L. In Fig. 
Etc) we show two T(E) curves for a wider electrode, keep- 
ing the length of the constriction fixed to L = 3a and 
changing the aspect ratio. The transmission correspond- 
ing to the lowest subband \E\ < Ai, is significantly re- 
duced with respect to the ideal ribbon T = 1, but remains 
finite (T(0) ~ 0.18) even for the narrowest and longest 
constriction. 

In order to highlight the peculiar properties of 
graphene, we give two arguments, which turn out to be 
inapplicable, to expect a vanishing transmission at zero 
energy. First, since low energy transmission is associ- 
ated to edge states and given that the edges of the elec- 
trode and the constriction are not connected, the trans- 
mission would vanish as iV c becomes much smaller than 
N. This argument fails because the current density of 
edge states has a sizable contibution in the center of the 
ribbon, which is smoothly connected to the constrictions 
considered in Fig. © . The second argument is based on 
the behavior of a square constriction with parabolic-band 
electrons. In that case, the different width of the semi- 
infinite ribbon W r and the constriction W c yields differ- 
ent band minima so that an energy gap is created in the 
narrow region. This can be modeled with square-shaped 

barrier of height V = ^ C^j ~ ^tt) and length L, m* 

being the effective mass of the electron. The zero energy 
transmission of a square barrier is always exponentially 
vanishing with the length, in contrast to our results. This 
second argument to expect a vanishing transmission at 
zero energy also fails because the situation for graphene 
ribbons is very different: Both wide and narrow sectors 
of the system have zero energy edge states. In other 
words, it is not possible to create a gap in the low en- 
ergy region by changing the width of the system so that 
it is not possible to create a barrier for the carriers in 
the lowest energy sub-band by geometrical means. In- 
terestingly, the transmission through a square potential 
has been calculated, in the case of infinitely wide rib- 
bon, and the transmission does not vanish for electrons 
incident perpendicular to the barrier—. 

The results in Fig. EJa) and (b) suggest that the low 
energy transmission is an increasing function of the con- 
striction width iV c and and fairly insensitive to constric- 
tion length L. In Fig. HJd) we plot the transmission at 
a fixed energy E = 0.05t as a function of N c /N for three 
families of constrictions. Square symbols correspond to 
a wide ribbon with N — 48, and L = 3a, whereas solid 
(open) circle symbols correspond to a narrow ribbon with 
N = 24 and L = a (L = 3a). It is apparent that as the 
aspect ratio increases, the transmission increases in av- 
erage, with super-imposed oscillations. This behavior is 
different from the analogous curve for transmission as a 
function of the barrier height T(E, Vq) for parabolic-band 
electrons in the tunneling regime (E < Vq), where the 
curve does not present oscillations. In the E > Vq regime 
the T(E, Vq) curve can present oscillations 24 . Therefore, 
the T(N C /N) curve looks closer to the latter case, even 
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FIG. 6: (Color On-line) Square constrictions, (a) T(E) 
curves for the N = 24, N c = 12 and two values of L. (b) 
T(E) curves for the N = 24, N c = 8 and two values of L. (c) 
T(E) curves for iV = 48, L = 3a for two constrictions with 
different 7V C . (d) T(E = 0.05i) as a function of the aspect 
ratio of the constriction N c /N for 3 families of constrictions. 



for very small energy. This is related to the lack of a gap 
in the density of states. It is also interesting that for both 
wide and narrow electrodes the transmission saturates to 
1 for the same value of N c /N ~ 0.7. As a final remark, 
only in the case of atomically narrow constrictions the 
transmission becomes almost zero. In summary, square 
constrictions are not very efficient in blocking electronic 
transport in zig-zag graphene ribbons. 



V. WEDGE-SHAPED CONSTRICTIONS 

We now turn our attention to a wedge-shaped constric- 
tion, shown in Fig. [7| As opposed to the square con- 
strictions considered above, wedge-shaped constrictions 
are always delimited by zig-zag edges. It can also be ex- 
pected that wedge-shaped constrictions are more stable 
than square constrictions, although further work should 
clarify this point. The structures considered have a top- 
bottom symmetry, in contrast to the asymmetric struc- 
tures considered in earlier work 21 * 2 ^. Naively, the wedge- 
shaped constriction would allow for adiabatic transport 
since the zig-zag edge gets never interrupted. In con- 
trast, we see that for ribbons with N = 48 atoms the 
constrictions shown in Fig. [7| yield a vanishing trans- 
mission in the whole low energy sector \E\ < A±. Wc 
have verified that these results are robust against the 
removal of a single atom from one of the edges. Surpris- 
ingly, these smooth nano-constrictions seem to be more 
disruptive for low energy transmission than the abrupt 
nano-constrictions considered in the previous section. 

Interestingly, the vanishing transmission is not associ- 
ated to a vanishing density of states in the device region, 
as shown in Figs. |H[c) and (d). The peaks in the DOS 
in the energy region with zero transmission is associated 



FIG. 7: Schematic atomic structure of wedge-shaped con- 
strictions 
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with low energy states very weakly coupled, if at all, to 
the electrodes. The width of the peaks in the DOS is 
limited by the numerical broadening used in the calcu- 
lation. In Fig. [5] we show a contour map of the DOS 
at E = and E = 0.04£, the energy at which the first 
peak of the device DOS is located for the structure shown 
on the right in Fig. [7| The maps reveal localization on 
the four edges of the constrictions, reflecting the fourfold 
symmetry of the device. Notice that the LDOS vanishes 
at the boundary with the electrodes, which explains the 
vanishing transmission. This LDOS could be observed in 
STM experiments^^. 

Whereas the probability density across one of the four 
edges has a bell shape for E = 0, it has a node for the 
finite energy case. The properties of these low energy 
non-conducting states in this wedge-shaped structures 
are very similar to those of zero-dimensional confined 
states. Confinement in semiconductor heterostructures 
is associated with the existence of an energy gap. The 
absence of such a gap in graphene makes it necessary 
to look for different strategies to confine electrons^iSSiSS. 
In this regard, nano-constrictions like the ones studied in 
this section can behave like quantum dots and might per- 
mit the study of Coulomb Blockade and Kondo physics in 
graphene structures. The physical origin of these bound 
symmetric wedge states (SWS) might be related to the 
formation of Kekule vortex structures discussed in the 
case of asymmetric wedge states£i*2£. In the asymmetric 
case, the suppression of the transmission occurs for very 
narrow energy windows, since one of the edges is not 
perturbed. It would be desirable to study whether the 
vanishing tranmission that we have found is related to 
the vanishing tranmission obtained analytically^! for the 
square barrier potential, in infinitely wide ribbons, for 
incidence angles different from zero. In summary, sym- 
metric wedge constrictions result in a gap in the trans- 
mission curve for \E\ < Ai, yet with a finite density of 
states in that interval, featuring very narrow peaks that 
mimic a discrete spectrum of confined states. 



FIG. 8: (Color On-line), (a) and (b) T(E) for the structures 
shown in Fig. [7| (blue) . (c) and (d) Corresponding density of 
states projected on the whole constriction region (red). 



VI. DISCUSSION AND CONCLUSIONS 

Here we discuss the validity of the approximations as- 
sumed for our calculations. Real samples could present 
several features absent in the idealized ribbon consid- 
ered here. On one side, the presence of chemical impu- 
rities, like hydrogen both in the edge and in bulk, water 
molecules, oxygen, etc., will affect the electronic struc- 
ture and transport properties of the system. In a suffi- 
ciently clean and ultra-high vacuum environment the ef- 
fect of chemical impurities could be negligible. In experi- 
ments, the graphene ribbon is deposited on an insulating 
substrate, which has not been included in our calcula- 
tion. First-principles calculations indicate that the effect 
of the substrate is weak 30 . 

The use of a single-particle model certainly fails, in 
a trivial sense, if the there are deviations from charge 
neutrality which would make it necessary to perform a 
self-consistent calculation including, at least, the Hartree 
contribution. This is certainly achieved by density func- 
tional calculations. We have verified that, for the struc- 
tures considered here, the electronic density does not de- 
viate significantly from charge neutrality. The single par- 
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FIG. 9: Local density of states for the right wedge structure 
shown in Fig. [7|at energies E = (a) and E = 0.04£ (b). 
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ticle approach also fails if the electron liquid happens to 
be different from a Fermi liquid, in which the low energy 
quasiparticles have the same quantum numbers than the 
free electrons. Such a scenario has been considered for 
two dimensional graphene 31 . A priori, this is a serious 
issue in a one dimensional system, where the Fermi liq- 
uid state is not stable with respect to electron-electron 
interaction and a Luttinger liquid is expected. The same 
statement applies for nanotubes. However, the single- 
particle Fermi liquid picture describes most of the ex- 
perimental results in nanotubes and the same can be ex- 
pected for graphene ribbons. 

We have used the same on-site energy (0) and hopping 
t for the edge and the inner atoms. This is an approxi- 
mation, but the results do not change significantly in the 
case of ideal ribbons**. Another limitation of our model 
is the one orbital approximation. Ab-initio calculations, 
including all the atomic orbitals, show that the dangling 
bonds present on the edge atoms form a flat band close 
to the Fermi energy^ 3 .. We have also ignored the spin de- 
gree of freedom. The flat band at zero energy is expected 
to spin split due to spontaneous magnetization induced 
by the Coulomb repulsior&S. This interesting issue de- 
serves more theoretical and experimental work. We have 
also ignored spin-orbit interaction, which is very small in 
carbon, although it has attracted some interest 3 ^ 3 ^, in 
part due to the Spin Hall effect. 

In summary, we have studied coherent transport in 
graphene nano-ribbons with zig-zag edges. The elec- 
tronic structure of the ribbons is described with a simple 
one orbital tight-binding approximation. We have fo- 
cused on narrow ribbons (W < 5nm ) for which only the 
lowest energy sub-band is expected to be experimentally 



relevant. Our results could be summarized as follows: i) 
The low-energy transport properties are robust with re- 
spect to isolated vacancies on the edges, ii) These are also 
robust with respect to square-shaped nano-constrictions. 
Linear conductance survives in most cases except for very 
narrow constrictions. This is at odds with the behav- 
ior of parabolic-band electrons in similar constrictions 
and reflects two non-trivial features of edge states: Their 
minimum energy is independent of the ribbon width and 
their current density profile spreads beyond the edges of 
the ribbon, in contrast with their density profile. Hi) 
In contrast to square constrictions, wedge-shaped nano- 
constrictions result in a gap in the transmission and re- 
sult in a zero linear conductance which is related to the 
appearance of localized low energy edge states. These 
edge states have properties similar to those of confined 
states in zero dimensions. Therefore, graphene wedge 
shape constrictions might have properties analogous to 
those of semiconductor quantum dots. Extensions of this 
work will address the spin degree of freedom, topological 
defects^ and armchair edges. 

We acknowledge useful discussions with M. Vozme- 
diano, A. Cortijo and B. Valenzuela. This work 
has been financially supported by MEC-Spain (Grants 
FIS200402356, MAT2005-07369-C03-01 and the Ramon 
y Cajal Program) and by Generalitat Valenciana (GV05- 
152). This work has been partly funded by FEDER 
funds. 

Note added: After the completion of this work we have 
been aware of a related theory paper— with transport 
calculations in wedge-shape nano-constrictions similar to 
that of figure J7J). 
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